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We investigate a four-impurity Anderson model where localized orbitals are located at vertices of a 
regular tetrahedron and find a novel fixed point in addition to the ordinary Fermi liquid phase. That 
is characterized by unscreened scalar chirality of a tetrahedron. In this phase, parafermions emerges 
in the excitation spectrum and quasiparticle mass diverges as l/|Tln T\ at low temperatures (T). 
The diverging effective mass is a manifestation of singular Fermi liquid states as in the underscreened 
Kondo problem. Between the two phases, our Monte Carlo results show the existence of a non 
Fermi liquid critical point where the Kondo effects and the intersite antiferromagnetic interactions 
are valanced. Singular behaviors are prominent in the dynamics and we find that the frequency 
dependence of the self-energy is the marginal Fermi liquid like — ImE ~ \uj\. 

PACS numbers: 72.15.Qm, 74.40.Kb, 75.10.Jm 
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I. INTRODUCTION 

The Kondo problem of multiple impurities has at- 
tracted much attention since the pioneering work of the 
two-impurity system about 25 years ago.— The presence 
of two competing energy scales, the Kondo temperature 
and the intersite exchange interaction, provides a nice 
play ground of quantum phase transition at zero temper- 
ature. Around the quantum critical point, non Fermi liq- 
uid (NFL) behaviors emerge and thermodynamic quan- 
tities and their dynamics exhibit singular temperature 
dependence. These physics have a connection to quan- 
tum criticality in heavy fermion systems^ and, thus, un- 
derstanding much simpler cluster problems is a key to 
capturing the nature of more complicated problems. 

Quantum fluctuations in such strongly correlated sys- 
tems can be controlled by tuning geometrical frustration. 
This is the case not only in spin systems but also in 
itinerant electron systems and the frustration leads to 
novel ground states such as spin liquid states As for 
the Kondo problem, frustration can be implemented by 
introducing a cluster of multiple impurities, and also con- 
trolling its geometry. When the cluster is frustrated, 
its ground state is degenerate due to nonmagnetic de- 
grees of freedom. In the simplest case where the clus- 
ter is a triangle^ the ground state is degenerate in the 
spin and the chirality sectors. Its Kondo problem was 
studied by the combination of the numerical renormal- 
ization group (NRG) and boundary conformal field the- 
ory (BCFT)^ and also by the continuous time quantum 
Monte Carlo (CTQMC) method^ and the renormaliza- 
tion group (RG)i^ It has been established that the frus- 
tration leads to a NFL phase. This phase is stable against 
perturbations that keep the triangular symmetry such as 
magnetic fields and particle-hole asymmetry^ which con- 
trasts with the unstable NFL in the two-impurity casei 

This paper reports extensive analysis on a minimal 
three-dimensional cluster problem with geometrical frus- 
tration, the Kondo problem of a regular tetrahedron. 



Tetrahedrons are simplex, i.e., the basic structure for 
three dimensional "frustrated" materials such as py- 
rochlore and spinel compounds. The result presented 
here would be helpful for further studies on the lattice 
problem in addition to the standpoint as a fundamental 
model analysis. 

This paper is organized as follows. In Sec. [TTl we intro- 
duce a microscopic model for four impurities. We analyze 
its two stable fixed points in Sec. IIII1 with emphasis on 
an unscreened phase. Section HVl is devoted to clarifying 
the global phase diagram and discussions about a criti- 
cal point between the two stable fixed points. Finally, we 
summarize the present results in Sec. |V] 



II. MODEL 

We first introduce a four-impurity Anderson model, in 
which localized s-wave orbitals are located at vertices of 
a regular tetrahedron and hybridize with a conduction 
band: 



H 



(1) 
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where tjj^ creates a conduction electron with wave vector 
k and spin a, d\ a creates a localized electron at position 
Xi(i = 1,2,3 or 4), n ?(T = d^d^ and other parameters 
are conventional ones. We also introduce the d-electron 
spin operator at Xj, for later use. 

It is useful to transform ip to partial waves^ that couple 
with d electrons. In the new basis, Hamiltonian ([1]) is 
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written as 



H 



J dk^ekipl„ a ipkao- + Vgjdl.^kaa- + h.c.) lo^S^t Kondo singlet 



Scalar chirality 

X=<SrSj*S, 




U terms, 



(2) 



where the orbital index a runs over four channels: 
(s, x, y, z). The orbital s(x, y, z) is the basis of the Ai(T 2 ) 
representation in the Td point group. The dependence 
of V a on wave number /c=|k| is ignored and replaced 
by the value at the Fermi wave number kp: V S =V[1 + 
3sin(fci?a)//cFa] 1 / 2 and V Xty , z =V[l — sin(fc f a)//^] 1 ^ 2 , 
with a=|xi — X2I being the impurity-impurity distance. 
V is a parameter proportional to Vq. The orbital energy 
is given as e s =e — 3t and e x ,y,z z 



-e + t. The new basis for 



d electrons is given by 

daa 77 ^ ] (^* )i^icrj 



(3) 



with F=(l, 1,1,1), P=(l, -1,1,-1), f»=(l, -1,-1,1), 
and f z =(l,l,— 1,— 1), and similar expressions for the 
conduction electrons. The bandwidth is set to 2D cen- 
tered at the Fermi level and the density of states is set 
to constant p = 1/2D for all ipkaa- 



III. TWO STABLE FIXED POINTS 

First, we investigate two stable fixed points in the 
Kondo regime, where the d electron is nearly half 
filling X)cr n icr~l- Note that there are two compet- 
ing energy scales. One is the on-site Kondo energy 
TK^Dexp(— |e|/4pVQ 2 ). The other is the inter-site ex- 
change interaction J, which includes superexchange and 
Ruderman-Kittcl-Kasuya-Yosida (RKKY) interactions. 
For Tk^-J, the Fermi liquid (FL) fixed point is stable, 
where each d electron forms a spin singlet with conduc- 
tion electrons. In the opposite limit, J3>T/<, four local- 
ized spins form spin singlets with a nonmagnetic double 
degeneracy: the E representation in the Td group. The 
two states in the doublet can be characterized as eigen- 
states of the scalar chirality (SC) 

X = Si ■ (Sj x S fc ), (4) 

with i j ^ k& The question is what happens when the 
interactions are switched on between the chirality dou- 
blet and conduction electrons? We show that conduction 
electrons cannot screen this doublet. Between the two 
stable fixed points, there exists a critical point as shown 
in Fig. [TJ and we discuss this in Sec. [IVJ 
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FIG. 1: (Color online) Schematic renormalization flow in the 
Kondo regime as a function of intersite exchange interaction 
J. 



represent it by a pseudospin 1/2, whose bases are 
1 



lit) 



-(i + /)(|UU)-|tUt», 



m = ^(i + /)[itm>-^(itm) 



(5) 

mt))]-(6) 



Here, \a\ oi<jz 04} represents a state with the z com- 
ponent of spin at at site i and / is the spin inver- 
sion operator. The effective model is strongly restricted 
by its symmetry. Since possible operators of the local 
E doublet are E®E=Ai®A 2 (BE, nontrivial couplings 
are those in the A 2 and E sectors. Taking into ac- 
count the classifications for the conduction electron bilin- 
ear operators, T 2 ®T 2 =A 1 @E®T 1 ®T 2l Ai®Ai=Ai, and 
A\®T 2 =T 2 , only a coupling in the E sector is possible: 



n int =g(T z Q + r x T), 



(7) 



where g cx i|l/i;| 2 /|e| 2 for £/— s-oo and (r z , t x ) are the spin- 
1/2 operators for the pseudospin with the E representa- 
tion. The conduction electron parts, 
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T 



V2 



(Afy -M X ), 



are the local orbital density with 



N a = ^2 j ' dkdk'ip\ aa ^ k > aa . 



(8) 
(9) 

(10) 



A. Unscreened phase 

1. An effective model 

In order to investigate the limit J Tk, we first derive 
an effective model of the local ground-state doublet. We 



They are the E representation and the "quadrupole" 
operator with [2,7"] = 0. This commutativity distin- 
guishes Eq. ([7]) from a conventional quadrupole Kondo 
model* and leads novel behaviors as discussed below. 
Note that Af s does not couple and the SC is represented 
by x~V^ Ty : which is the A 2 representation and cannot 
appear in the Hamiltonian. 
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2. Perturbative renormalization group 

Next, we analyze the effective Hamiltonian ((TJ) by RG. 
It is remarkable that the one-loop order term vanishes 
and the leading order is the two-loop level, and the RG 
equation is 



(ii) 



The absence of 0{g 2 ) terms in the RG equation is traced 
in [Q,T] = 0. Solving the RG equation, we obtain 
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P9cs 



^\n(T*/D cS ) 



(12) 



with T*=D exp [1 / (pg) 2 ] as the effective bandwidth is de- 
creased to D e g. This asymptotic logarithmic form is sim- 
ilar to that in the underscreened Kondo models aside 
from the extra power 1/2. Thus, the coupling constant 
approaches zero with lowering temperature (energy). 
The impurity degrees of freedom cannot be screened and 
there remains a residual entropy In 2. The pseudospin 
susceptibility has the Curie like temperature dependence 



1/T 



(T T T*(r)T*(0))dT~T- 1 , (13) 



where T T represents time-ordering. Summing up the 
leading logarithmic terms, we find that the impurity spe- 
cific heat coefficient Ci mp /T diverges as 



1 



T Tin (T*/T) 



(14) 



This means a diverging effective mass, and the power 3 
is different from that in "singular" Fermi liquid states.— 



3. Boundary conformal field theory and numerical 
renormalization group analysis 

The presence of the unscreened phase is explicitly 
demonstrated in our NRG^ 1 - calculations. Energy lev- 
els have degeneracy expected for their quantum numbers 
(see below) plus an additional factor 2. This double de- 
generacy is a manifestation of the unscreened SC. The 
quantum numbers of each level consist of three parts: 
total charge Q, spin j, and the orbital part. The en- 
ergies and the quantum numbers arc in complete agree- 
ment with the prediction below derived via "phase shift" 
in the orbital sector in the BCFT language^ assuming 
an unscreened SC. This confirms that the unscreened lo- 
cal object in this phase is a SC, or equivalently different 
configurations of two spin-singlet pairs£ 

Let us explain our BCFT analysis. We first in- 
troduce a conformal embedding suitable to Eq. (0. 
The free Hamiltonian is represented by in addition to 
the global U(l) charge current, the spin and orbital 



currents following the Kac-Moody algebra SU(2)3 and 
SU(3)2, respectively. One possible conformal embedding 
is U(l)<g>SU(2)3<g>SU(3)2, and the energy eigenvalues of 
the free Hamiltonian are represented as 
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i(i + i) . c su(3) 



integer 



, (15) 



where we have defined the effective onc-dimcnsional space 
in [—1,1], vf is the Fermi velocity, and csu(3) is the 
eigenvalue of the Casimir operator in the SU(3) sector: 
csu(3) = 0_for 1, c su(3) = 4/3 for 3 and 3, c su(3) = 10/3 
for 6 and 6, and Cgu(3) = 3 for 8. 

The interaction ([7]) preserves each of the charges in the 
T2 orbital separately, although it breaks the orbital SU(3) 
symmetry. This means there remains two conserved U(l) 
charges in the orbital sector: 



Q 
T 



1 
1 

72 



(J\fy -ST X ), 



with 



(16) 
(17) 

(18) 



Note that the total orbital charge J\[ a differs from the 
local density N a and Q = J2 a N a - Factorized out with 
these two U(l) charges, the CFT for the remaining or- 
bital part SU(3) 2 /[U(1)] 2 is known as the parafermion 
(PF) one with the central charge c = 6/5.— Thus, the 
energy spectra for the free Hamiltonian are written as 



E = — 



91 
12 



jU + 1) 



h A + integer 

4 



(19) 



where TZ 2 = T 2 + Q 2 and A is the dimension of the 
primary fields in the PF sector (see the caption in Tabic 
UJ. Equations (|15|) and (Tl"9l give exactly the same spectra. 
Note that there is 0(2) symmetry between Q and T. Q 
and TZ terms in Eq. (|T9l) can be also represented as 
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H 2 _ M 2 +M 2 +M 2 



(20) 



which is symmetric with respect to the orbital indices, as 
it should be. 

Tabic U and Figf^a) show NRG spectra obtained with 
keeping 10 4 states at each step. We have checked that the 
results do not change when 10 4 more states are added. 
We find that the energy eigenvalues Em at the Nth. RG 
step are given by the free spectrum modified by an addi- 
tional potential scattering in the 0(2) sector: 



I n 

E N oc E° ± g N K, 

ttvf 



(21) 
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TABLE I: Quantum numbers, energy, and degeneracy of 128 
low-energy states for odd N. In the parafermion (PF) sector, 
following notations in Ref.— , there are eight primary fields: 
{1, af, <7|_, 03, ipi,ij)2, 1P12, p}- cr's and ijj's are abbreviated sim- 
ply as <j and ip, respectively. Their scaling dimension A is 
indicated in parentheses. The energy is measured from the 
value of the non-interacting ground states and then scaled 
appropriately. 
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FIG. 2: (color online), (a) NRG spectra for odd N for g = 
0.5D and A = 3. At each N, the energy of the lowest-energy 
states are subtracted. States are labeled by Tl 2 and the sign 
of the energy listed in Table [I] (b) En as a function of 



space as 

C = V £n = i J2 - [QnT-n - T n Q- n ] + C . (23) 

n>0 n>0 

Here, Q n and T n are the nth Fourier modes for the U(l) 
bosons and 



£0 — Qo&t — Th^Qi 



(24) 



with <7tv constant. Note that this form remains 0(2) 
symmetric and is derived from the free spectra by shifting 
K -> TZ±2g N in Eq. ((TSJ), i.e., just a "phase shift." The 
factor ± indicates unscreened SC as in the Ising Kondo 
caseJ^ We note that the simple "phase shift" leads to PF 
excitation spectra for g^ > 0. 

The Nth step in the NRG calculation is related to the 
energy scale DA~ N / 2 with discretization parameter— A 
and we set A = 3. Replacing this for D c g in the two- 
loop expression of g c s, we find g^ oc 1/yN. Indeed, we 
obtain 



En — Eqo oc ± 



1 



N 



n. 



(22) 



in the NRG results as shown in Fig(5] (b), which con- 
firms the two-loop results. This 1/\N dependence is a 
clear contrast to 1/N dependence in underscreened and 
ferromagnetic Kondo models. 



4- A hidden angular momentum 

We have observed that the 0(2) symmetry in the or- 
bital space Q-T is preserved in the energy spectrum of 
the interacting system g ^ 0. This holds asymptotically 
where the bosonization becomes exact in the low-energy 
limit. We demonstrate this below by explicit calcula- 
tions. 

The 0(2) symmetry is related to a hidden conserved 
quantity of angular momentum. This angular momen- 
tum is nontrivial and defined in the two dimensional Q-T 



where 0q(0t) is a canonical conjugate field to Qo(7o), 
with [Q , 9q] = [Tq, 9f]=i. C is conserved for g = 0, but 
is not for g ^ 0. The conserved one for g 7^ is 



J = C 



1 



(25) 



which satisfies ["Hint, J] = 0. Wc find that eigenvalues of 
C are all integers, since S y — Xm>o w i* n = 
constituting the SU(2) algebra together with the partners 



S l — E n >0 
dz 

°n 

for n > 0, and 



x, z), where S^'s are given as 



in 
1 

2ri 



[QnQ-n - T n T- n + (n 
(Q n T- n +T n Q-n), 



s z = - A {Ql-T 2 + el-e 2 r ), 

S* = \(QoT + 9 Q e T ). 



-n)}, (26) 
(27) 

(28) 
(29) 



S^s satisfy S£j = ie ljk S nm S k , and this leads to the 
SU(2) algebra of [S\ &} = ie^ k S k . Thus, the eigenvalues 
of C are all integer, and C is, indeed, the "orbital" angular 
momentum. From this, the eigenvalues of J should be 
all half integers. Noticing that the total Hamiltonian 
is invariant under J^r —J, e.g., f|"<H-4 and Q-^—Q, all 
the energy eigenstates should be degenerate, which is the 
direct consequence of the unscreened state. 
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B. Fermi liquid phase 

For Tj{ 3> J, the situation is conventional and the 
ground state is a FL. Each ci-electron spin forms a spin 
singlet with conduction electrons. This phase is contin- 
uously connected to the ground state of the spin 5 = 2 
four-channel Kondo problem. The way of fully screened 
processes are similar to those in a Kondo singlet phase in 
the two-impurity Kondo model, where it is essentially an 
5 = 1 two-channel modeli Since there is an asymmetry 
between A\ and T2 orbitals in Eq. ([2]), screening takes 
place in two stages as T decreases. Indeed, we have con- 
firmed the Fermi liquid properties by using Monte Carlo 
simulation as is explained in the next section. 

IV. CONTINUOUS-TIME QUANTUM MONTE 
CARLO ANALYSIS 

Now, we unveil the NFL properties and determine the 
global phase diagram of this model by using CTQMC— 
All the data presented below are for U=—2e=1.5D, 
t=Q.2D, and &i?a=3. With these values, J is dominated 
by the superexchange ~4t 2 /t/ with a small antiferromag- 
netic contribution of the RKKY term. For other values 
of k,Fd, the results are qualitatively the same as long as 
J is antiferromagnetic. 

A. Thermodynamics 

FigureEJa) shows the susceptibility of local pseudospin 
X zz for several values of V. For V < 0.275D, the chirality 
exhibits unscreened behaviors Xzz^T _1 at low T . which 
is consistent with the discussions above. For larger V, 
Tk 3> J, and thus the ground state is a FL. There, %zz is 
small and constant at low T. Figure GO^b) shows that the 
double occupancy (nifHij,) ~ 1/2 — 2(5* ) increases as V 
increases. This implies that the d electrons become more 
itinerant as V increases. The intersite spin correlation 
(5*5f) is large and antiferromagnetic for small V, while 
it is suppressed for large V. As for local configurations 
in the four-electron sector, the occupancy of the ground- 
state doublet E with spin 5 = decreases with increasing 
V, while those for the first and the second excited states 
T2(5 = 1) and Ai(S = 2) increase as shown in FigJSJc). 
This is a natural consequence of the Kondo screening, 
which mixes local configurations with different 5. 

B. Single-electron dynamics 

A drastic change also appears in the single-electron dy- 
namics upon varying V. Figure0ta) shows the imaginary 
part of the electron self-energy H(iui n ) for the impurity 
T2 orbital, for T = 0.005-D as a function of the Matsub- 
ara frequency tu n . For large V > V* ~ 0.3D, lmT,(iuj n ) 
is linear in w n for small uj n , as expected for the FL state. 




0.2 0.4 0.6 



V/D 




V/D 



FIG. 3: (Color online), (a) T dependence of susceptibilities 
Xz Z for V/D = 0.2 - 0.4, U = -2e = 1.5D, and t = 0.2D. 
(b) Equal-time correlation functions: {Si 2 }, (nifnij,), and 
(Si Sf). ( c ) V dependence of the occupancies of the E, T2 and 
Ai states (per spin and orbital) in the four-electron sector. 



For smaller V, it shows a diverging behavior instead and 
this indicates that the T2 electrons are localized in the un- 
screened SC state. This is consistent with our conclusion 
based on the exchange model (J7J), since the localization of 
the T2 electron leads to the absence of resonance peak at 
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FIG. 4: (Color online). ui n dependences of (a) lmE(iuj„) and 
(b) ReE(iaj n ) for T = 0.005 D. Line s for V > 0.3D in panel 
(b) represent fit by ci + C2V8' 2 + oj% and the values of Ci, C2, 
and S are shown in the inset. 



V > V* is well fitted by - a + c 2y / S 2 +u% with 5>0. 
As shown in the inset, the critical point is realized by 
<5 — s- 0. This w„ dependence leads to the NFL form of the 
retarded self-energy: 



ImS R (w) 



(30) 



which is the same as in the marginal FL theory.— The 
NFL point V* ~ 0.3D is consistent with the transition 
point of V in FigH Thus, our CTQMC results for fi- 
nite temperatures strongly suggest that there is a quan- 
tum critical point separating the FL and the unscreened 
phases at T = 0. 

As for the Ai electron, the imaginary part of the self- 
energy shows FL behavior for all the values of V we ex- 
amined. This persists for smaller t and is related to the 
fact that the A\ orbital cannot interact with the non- 
magnetic E doublet apart from potential scattering. 

This NFL is stable against various perturbations which 
are relevant at some other types of critical points. For ex- 
ample, it is stable against magnetic field and particle-hole 
asymmetry, and this is understood from the difference in 
the ground-state entropy in the two stable phases. Dis- 
tortion that breaks the symmetry may give nontrivial 
effects on the stability of the NFL. Examination of the 
stability of the NFL against distortions and clarifying the 
NFL finite-size spectra are our future work. 



V. SUMMARY 

We have investigated an Anderson model with four 
impurities on a regular tetrahedron. We have found that 
the system has two stable fixed points. One is the new 
fixed point where an emerged SC of the tetrahedron is not 
screened by conduction electrons for small V with diverg- 
ing effective mass and the emergent PF excitations, and 
is a new class of singular FL states. The other is a Kondo 
screened FL state for large V. Our CTQMC results have 
revealed that the new critical NFL state appears, ac- 
companying the d-electron localization transition at zero 
temperature and the self-energy exhibits marginal FL- 
likc frequency dependence. 



the Fermi level in the conduction-electron Green's func- 
tion for the exchange model ([7]). The crossover around 
V* is similar to that near a critical end point of metal- 
insulator transitions in strongly correlated systems, e.g., 
Mott's transition^ 

The real part RcE(ia; rl ) shows an even more peculiar 
behavior near V* as shown in FigQJb). ReE(io; rl ) for 
V ~ V* shows a nonanalytic behavior, ~ c\ + C2|w n |, 
where c\ and C2 are constants. The w„ dependence for 
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